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Abstract 

O . We study the effective action in Euclidean Yang-Mills theory with 

a compact simple gauge group in one-loop approximation assuming a 
covariantly constant gauge field strength as a background. For groups 
^ . of higher rank and spacetimes of higher dimensions such field config- 

urations have many independent color components taking values in 
Cartan subalgebra and many "magnetic fields" in each color compo- 
nent. In our previous investigation it was shown that such background 
is stable in dimensions higher than four provided the amplitudes of 
"magnetic fields" do not differ much from each other. In present paper 
we calculate exactly the relevant zeta-functions in the case of equal 
amplitudes of "magnetic fields" . For the case of two "magnetic fields" 
with equal amplitudes the behavior of the effective action is studied 
in detail. It is shown that in dimensions d = 4,5,6,7 (mod 8), the 
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perturbative vacuum is metastable, i.e., it is stable in perturbation 
theory but the effective action is not bounded from below, whereas 
in dimensions d = 9, 10, 11 (mod 8) the perturbative vacuum is abso- 
lutely stable. In dimensions d = 8 (mod 8) the perturbative vacuum is 
stable for small values of coupling constant but becomes unstable for 
large coupling constant leading to the formation of a non-perturbative 
stable vacuum with nonvanishing "magnetic fields" . The critical value 
of the coupling constant and the amplitudes of the vacuum "magnetic 
fields" are evaluated exactly. 



PACS numbers: ll.lOKk, 11.15Tk, 11.15.-q, 12.38Aw, 12.38Lg 
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1 Introduction 

The low-energy behavior of the non-Abelian gauge field theories is an open 
problem of quantum field theory. At low energies the effective interactions 
become strong and break down the usual perturbation theory, which works 
very well in the high-energy limit. This means that the vacuum at low 
energies has a far more complicated structure than the perturbative high- 
energy one. 

To analyze the infra-red behavior of a SU(2) Yang-Mills theory Savvidy 
[p]] put forward an explicit ansatz for the vacuum gauge fields in form of a 
covariantly constant gauge field strength. Such a field configuration is neces- 
sarily Abelian, i.e., it takes values in the Cartan subalgebra, and, hence, since 
rankS'L r (2) = 1, has the only nonvanishing color component. On the other 
hand this component has the form of a constant magnetic field. Therefore 
such a configuration is usually called constant chromomagnetic field. Savvidy 
showed that due to the quantum fluctuations of the gauge fields the energy of 
such a field configuration is below the perturbative vacuum level, which leads 
to infra-red instability of the perturbative vacuum under creation of the con- 
stant chromomagnetic field. Further investigations @, [3|, |, [5|, ||, [7], || showed 
that the constant chromomagnetic vacuum itself is unstable too, meaning 
that the real physical nonperturbative vacuum has even more complicated 
structure. 

In our recent papers || [UJ we extended the Sawidy's investigation by 
considering more complicated gauge groups and spacetimes of dimensions 
higher than four. We showed that there exist more general nontrivial field 
configurations that turn out to be stable. 

In present paper we are going to consider some explicit example of such 
background field configurations when the problem can be solved exactly to 
the end. 

2 Infra-red regularization of the low-energy 
effective action 

We consider the pure Yang-Mills model with an arbitrary simple compact 
gauge group in an Euclidean (/-dimensional flat spacetime. The classical 
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action of the model is 

S = J dx^tr M F^, (1) 

where g is the coupling constant, T^ v is the gauge field strength taking values 
in the Lie algebra of the gauge group and tr^ means the trace in adjoint 
representation. 

In the case of covariantly constant background fields, 

= 0, (2) 

the Euclidean effective action does not depend on the gauge || and is usually 
given in one-loop approximation by (see e.g. [O) 

T = S + hT (1) + 0(h 2 ), (3) 



r (1) = -logDet--logDet-^. (4) 

Here 'Det' means the functional determinant, fi is a dimensionful renormal- 
ization parameter and A and D are elliptic second-order differential operators 
of Laplace type 

A* = - 2^, (5) 

D = -□, (6) 

where □ = V^V 1 is the covariant Laplacian, the covariant derivative V M 
being taken in adjoint representation. 

In practice, to study the infra-red behavior of the system one has to 
be a bit more careful. Namely, to be able to carry out the intermediate 
calculations one should impose some infra-red regularization by introducing, 
say, a regularization mass parameter M and take it off at the very end. The 
regularization parameter M must be sufficiently large to ensure the infra- 
red convergence. Then one has to do the analytical continuation to the 
neighborhood of the point M = and to take the limit M — > 0. If there are 
no infra-red divergences then this procedure is harmless. In the case of non- 
trivial low-energy behavior there appears an imaginary part of the effective 
action or some infra-red logarithmic singularities, so that in addition to the 
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ultra-violet renormalization it could also be necessary to do the infra-red 
renormalization. 

The infra-red regularized effective action T(M) is not analytic, in general, 
at the vicinity of the point M — 0. Therefore, one should assume the infra- 
red regularization parameter to have an infinitesimal negative imaginary part, 
i.e., M 2 — > M 2 — ie. This enables one to choose the correct way of analytical 
continuation leading to the correct sign of the imaginary part of the effective 
action. Thus, strictly speaking, the effective action should be defined by 

r d ^ f lim.r(M), (7) 

M— *0— ie 

In the one-loop approximation one can define the infra-red regularized 
effective action by 

r (1)W i logDet (^!)_ logDet (^!), (8) 

and suppose M 2 to be large enough, so that the operators A + M 2 and 
D + M 2 are positive, i.e., do not have negative modes. 

The functional determinants are also known to be ultra-violet divergent. 
To regularize the ultra-violet divergences one can take any usual regulariza- 
tion, e.g. the dimensional regularization or the zeta-function one etc. One 
should stress that the meaning of the infra-red regularization is very differ- 
ent from the ultraviolet one. Although the ultraviolet regularization is purely 
formal method to get physically predictable results, the infra-red regulariza- 
tion parameter can take, in principle, a direct physical meaning, something 
like Aqcd- Anyway, the infra-red regularization should be taken off at the 
very end, after the ultraviolet regularization. 



3 Zeta— function and heat kernel 



Within the zeta-function regularization [12, 13, 14] the effective action has 
the form 



r ( i)(M) = --c YM (o,M), 



(9) 



where 
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Cym(p, M) = Ca+m 2 (p) - 2Cd+m2(p) 



(10) 



is the infra-red regularized total Yang-Mills zeta-f unction. 

The zeta-function of an elliptic differential operator L can be defined in 



terms of the associated heat kernel as follows [112;, O, WA 



2p 00 

CM = H 2p TrL~P = Jdx^p - J dtlThTU L (t) 



r(p) u 

where 'Tr' means the functional trace, 'tr' is the usual matrix trace and 



(11) 



U L (t) = ejqp(-tL)5(x,x') 



is the heat kernel diagonal. Therefore, the Yang-Mills zeta-function is ex- 
pressed in terms of the heat kernel diagonals of the operators A and D 



where 



Cym(p,M) = Jdx^J dtt^e-^Uyuit), 



(13) 



U YU {t) = tr Ad {tr r U A (t) - 2U D (t)} 



(14) 



and 'try' means the trace over the vector indices. 

The heat kernel diagonals of Laplace type operators L = — □ +P on 
covariantly constant background are calculated in |§ using an algebraic ap- 
proach developed in [fT5| . This approach was developed further, also to the 
case of curved manifolds, in [16|, [L7|, |18| . (For review see |1| ^(| ^TJ). 

By using the results of M we have 



U YM (t) = (4vrt)- d / 2 tr Ad { det T ( — (tr T exp(2^) - 2 ) \, (15) 

where T = {J 7 ^} is a matrix with vector indices and 'defy' means the 
determinant over the vector indices. 

To calculate the trace over the group indices we observe first that covari- 
antly constant Yang-Mills fields are necessarily Abelian, 
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[J>,^a/3]=0, (16) 

and, therefore, take their values in the Cartan subalgebra. Thus the maximal 
number of independent fields is equal to the dimension of the Cartan subal- 
gebra, i.e., to the rank of the group r. Using the property (|I6| ) one obtains 



U YM (t) = (4nt)- d / 2 ^r(d-2) 
where 



2 det T ( - 



tF 



1/2 



,sin(£F), 
F = {F a >J, 



tr T cos(2tF) -2 



(17) 
(18) 



a a are the positive roots of the Lie algebra of the gauge group and the sum 
runs over all positive roots. The number of positive roots is p = (n — r)/2, 
where n is the dimension of the gauge group. 

Further, by orthogonal transformations one can always put the antisym- 
metric d x d matrix F in the block-diagonal form with q two-dimensional 
antisymmetric blocks on the diagonal, all others entries being zero 
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Visually one can think of such field configuration as a number of 'magnetic 
fields' in different directions and call the invariants Hi the amplitudes of 
'magnetic fields'. They can be expressed in terms of the basic invariants of 
the matrix F, 



by solving the equations 



+ 



< : 



H 2 [d /2] + 



-(-l) fc tr T F 



+ H l 



+ ff2[d/2] 



2k 



(20) 



h 



(21) 



l[d/2] 



One should note that, although the matrix F depends linearly on the 
roots, its invariants ffj depend on the roots in a very nontrivial way. 

By making use of flT9| ) one can compute the trace and the determinant 
over the vector indices in eq. fll?D and obtain finally 



U YM (t) = (4nt)- d ^^r(d-2) 
+2 S ) 



Hence the total zeta-function for the gauge fields, (|13D , equals 



rf-2 + 4 sinh2 ( tH i 

i<i<? 



" 2 5iSA sinh (^ 



(22) 



r u 2p °r ( 

( YM (P,M) = j dx(An)-^ 2 ^ I dtt^ d / 2 ~ 1 e~ tM2 \r{d-2) 

-) + sinh 2 (^)j|.(23) 



For Rep > d/2 the integral ( p3| ) over t converges at t — > 0. For sufficiently 
large M it also converges at t — > oo for arbitrary values of amplitudes of 
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'magnetic fields' Hi. It is easy to see, however, that for M = there are 
some field configurations that break down the convergence of this integral at 
t — > oo. The simplest example is the case of only one 'magnetic field', i.e., 

9=1- 

In || it was found a criterion of infra-red stability, i.e., the convergence of 
the integral ( p3|) over t at t — > oo in the limit M — > 0. Roughly speaking, the 
integral fl23|) converges when the amplitudes of 'magnetic fields' Hi do not 
differ much from each other. More precisely, the vacuum is infra-red stable, 
i.e., the integral ( p3|) converges at t — > oo for M = 0, when the background 
'magnetic fields' satisfy the condition of stability, 

i 

max < Hi for any root a, (24) 

!<*<<? ,7^ 
l<i<q 

where the prime at the sum means that the summation does not include the 
maximal invariant. 

Obviously, the condition (|24]) can be fulfilled only in the case when the 
number of 'magnetic fields' is not less than two, q > 2, and, hence, only in 
dimensions not less than four, d > 4. The case q = 2 in d = 4 is the critical 
one (see the detailed discussion below). That is why we also study in present 
paper the higher-dimensional case (d > 5) in detail. Even in this case the 
vacuum can be unstable. Since the stability is supported by approximately 
equal amplitudes of 'magnetic fields', we study in this paper the 'most stable' 
configuration when all the 'magnetic fields' have equal amplitudes: 

Hi = H, (i = l,...,g). (25) 

This assumes certain algebraic relations between the basic invariants Ik of 
the matrix F 

f j v 1/2 , x 1/4 ,j \ l/(2[d/2]) 

^K) =(f) =••■= W ■ (26) 

In this case the Yang-Mills zeta-function takes the form 



r ii 2 p °r ( 

(ym(p,M) = dxU7i)- d / 2 -^— dtt p - d / 2 - 1 e- tM2 \r(d-2) 

J r(p)y { 
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For sufficiently large M the integral (|27|) converges in the region Rep > 
d/2 and defines an analytic function. Therefore, changing the integration 
variable t — > t/H we get 



(ym(p,M) = J dx(4n)- d /*^ TiP T ^ /2) lr(d - 2)M d ~^ 

(28) 

where 

J Sjq (z) = (d-2)b s , q (z)+Aqs(s-l)b s - 2 , q -2(z), (29) 

i 00 -tz 

J «W = fH/* r " ,+ '^' (30) 

v ' 

Thus the regularized zeta-function is expressed in terms of the functions 
b s>q (z). These functions are studied in detail in the Appendix. It is well 
known that the coefficient at (—log fi) in the effective action is determined 
by the zeta-function at zero. From eq. fl28|) we immediately get 

Cym(0,M)=0 for odd d, (31) 

and 



Cym (0,M) = / d x(4.)^ f ^ y { r ( d -2)M^ + 2|: J d/2 „ ) H^} 

(32) 

for even d. 

Let us take the limit M — > — is here. The contribution of the 'free' 
fields is proportional to M d and vanishes in this limit. Defining 

J Siq = f lim J s , q (z) (33) 
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and using the formulas of the Appendix we find that J S)9 for odd positive 
s = 2k + 1 vanishes, 



J2k+i, q = 0, {k = 0,1,2,...)- 
whereas for even positive s = 2k 

J2k, q = {d- 2)a 2 k, q + 8qk(2k - l)a 2 k-2, q -2 
with (see Appendix) 



d\ 2k ( t 



a 2k, 



t=0 



and 



q \dt j Vsinht, 
Therefore, we have from ([32]) 

C YM (0, 0) = for odd d > 3 and odd d/2 > 1, 

2 



(34) 
(35) 

(36) 
(37) 



Cym(0,0)= / rfx^(4vr) 



-d/2 



T(d/2 



-?J d /2, q H d/2 for even d/2 > 2. (38) 



4 One-loop effective action 

Using the zeta-function ( pq) we obtain from (§) then the effective action: 



r (1) (M) = | 



'M 2 



a>0 



(39) 



for odd d 



and 



r (1) (M) 



2 V ; r(d/2 + l) 

M 2 



x { r (d - 2)M d log — - *(d/2 + 1) - C 



+2 x; # d/2 



o>0 



+ (log^-^d/2 + 1) 



- C 



(40) 
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for even <i, 



where 



Jd/2,q{ Z ) 



d J s,q{z) 



ds 



(41) 



=d/2 



z) = dlogT(z)/d z and C = — ^(1) is the Euler constant. 
Very similar formulas for the zeta-function and the effective action were 
obtained in general case in [p| in a bit different context (see especially the 
eqs. (2.25), (2.28) and (2.29) in Sect. 2 therein). 

Thus the effective action is determined by two quantities — Jd/2,q{M 2 / H) 
and J'd/2,q(M 2 1 H), i.e., by the values of the function J s , q (z) at the positive 
integer and half-integer points, s = k/2, and its derivative at integer points, 
s = k. Using the eq. (|29|) and the properties of the function b S) q(z) studied 
in the Appendix, it is not difficult to obtain all these coefficients. 

By taking off the infra-red regularization and using the eq. fl34|) we obtain 
finally 



(i) 



r tt(—\ V^- 1 )/ 2 

J dx E W"'^ p (d/ 2 + i) J */^ Hd/2 fOT odd d > 3 ' 



(42) 



and 



dx (47r) 

a>0 



-d/2 



r(d/2 + i) 



J'd/2, q H d/2 for odd d/2 > 1, (43) 



(i) 



^ ' r(d/2 + i) 



x 



H d/2 f 



log4-*W2 + l)-C 



(44) 



for even d/2 > 2. 

Thus one has to calculate the values of the function J Si 9 and its derivative 
J' s at half- integer points and integer points. We distinguish three essentially 
different cases: i) q — 1, ii) q — 2, and iii) g > 3. 
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4.1 One 'magnetic' field 

The case q — 1 is realizable in any dimensions d > 2. Therefore, we need to 
study the function J Sjl in the region Res > 2. It should be noted, however, 
that this is the only possible case in four- dimensional spacetime of Lorentzian 
signature 0. It is this case that was studied by Savvidy [0. 

Using the formulas of the Appendix we obtain from eqs. (^) and (|33|) in 
this case 

J s l = 4(d-2)(2 s - 1 -l)(27r)^cos (^s) T{s + l)C(s) + 2s (e~™ s + l) . (45) 
Therefrom we obtain for odd dimensions, i.e., for half-integer s = d/2, 



Jd/2,1 = (-l) [(d+1)/4l 2(d-2)(2 d/2 - 1 -l)v / 2(27r)- d / 2 r(rf/2 + l)C(c//2) 

+d-t(-l) {d - 1)/2 d, (46) 



for odd d > 3, i.e., d = 3,5,7, ... . 



For even dimensions, i.e., integer s = d/2, we find from (f45D in accordance 
with eqs. (|§ and 



J d/21 = for odd d/2 > l,i.e., d = 2,6,10,... 
and, by using |23[ 



((2k) = (-l) fc+ W fc ^fyy, 
where Bk are the Bernulli numbers, 

Jd/2,1 = ~2{d - 2){2 d ' 2 - 1 - l)B d/2 + 2d, 
for even d/2 > 2, i.e., d = 4, 8, 12, . . . . 
The derivative J' d ^ 2 1 for even dimensions is 



(47) 
(48) 

(49) 



j' 

J d/2, 1 



i-lY d ^(d - 2){2 d ^ - l)T(d/2 + 1)tJP^ + ™* (50) 
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for odd d/2 > 1, i.e., d = 2, 6, 10, . . . 



J d/2, 1 



-2{d - 2){2 d ' 2 ~ 1 - l)B d/2 
2<i/2— 1 

log2-log(27r) 



C'(d/2) 



2 d/2-l _ ! 



C(d/2) 
+ 4 - wrd 



+ *(d/2 + l) 



(51) 



for even d/2 > 2, (d = 4, 8, 12, . . .). 
Substituting (g§), (§§) and (0) in (g^-Q we find that the effective 



action gets a negative imaginary part 
Im T(x) = — J dx (4 



d/2 



for 9 = 1, 



(52) 



r(d/2)" 

which indicates on the instability of the chromomagnetic vacuum with one 
'magnetic field', i.e., q = 1. It is remarkable that this form is valid for any 
dimension. 



4.2 Two 'magnetic' fields 

Now let us consider the most interesting case of two 'magnetic' fields, i.e., 
q = 2, that is possible in (d > 4). Again we have to study the region 
Res > d/2 > 2. This case is critical because for q = 1 there is a strong 
infra-red instability, whilst for q > 3 the model is infra-red stable. The limit 
z — ► — ie of the function J Si i(z) is not regular. From eqs. (^9|) and ( j33|) 
and the Appendix we have now 



2 s+1 (d -2){s- l)(2n)- s cos ( ^-s ) F(s + l)C(s 



z— >0— ie 



71 



+8s(s - l)z 



s-2 



(53) 



Note that the limit z —>■ — ie is well defined only in the region Re s > 2, 
i.e., for d > 5. In this region we have 



J,, 2 = 2* +1 (d - 2)(s - l)(2vr)- s cos ( -s ) T(s + l)C(s). (54) 



7T 
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Therefore, 



Jd/2,2 = {-if^l^-^id - 2) 2 (2vr)- d / 2 r(rf/2 + l)C(<*/2), (55) 
for odd d > 5, (d = 5, 7, 9 . . .) 



and (in accordance with ([34]) and (j35l) ) 

2 = for odd d/2 > 3, (d = 6, 10, 14, . . .) (56) 

■V.2 = -2 d/2 - 1 (d - 2) 2 5 d/2 for even d/2 > 4, (d = 8, 12, 16, . . .)• (57) 
The derivative J' d , 2 2 for even dimensions is given by 

J' d/ 2,2 = (-l)^+ 2 )/ 4 ^/ 2 - 2 (d - 2) 2 r(d/2 + 1) (58) 
for odd d/2 > 3, (d = 6, 10, 14, . . .) 

and 



4/2,2 = -2 d/2 1 (d-2) 2 B d/2 



s ^ 1 f// 2 ' + (d/2 + 1) + log 2 - log(27r) + 2 



Cd/2 v ' ' m ' d-2 

(59) 

for even d/2 > 4, (d = 8, 12, 16, . . .)• 

The case d = 4 is not regular one because the point s = 2 is singular in 
the limit z — > 0. The value J 2j 2 is still finite 

4 44 

J 2 , 2 = -- + 16 = y, ford = 4, (60) 

where the second contribution, +16, comes from the z-dependent term in 
d53|). Note that this cannot be obtained from (|57D by putting d = 4, but is 
in fully accordance with eq. ([35|) for g = 2 and d = 4. 

However, the derivative J' d / 2 2 (z) for d = 4 diverges in the limit z — > 

4 2 (^)|_ = 16 log z - \ ( ^ - C - log(27r) + log 2 + 1 ] + 22 + o(z). 



3 VC(2) 



(61) 
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Therefore, from eq. ( fjQ|) we find that the effective action T(M) has a loga- 
rithmic infra-red divergency in the limit M — > 

r {1) (M)| M ^ o = -|rfxE(4vr)-Wlog^ + 0(l), (q = 2, d = 4). 

(62) 

4.3 More than two 'magnetic' fields 

The case g > 3 is realizable only in dimensions d > 6. Therefore, we need, in 
fact, to know J s ? only for Re s > 3. In this case J S}q (z) is an entire function 
of s for Re z > (2 — q). This means that it also remains entire function of s 
in the limit z — > 0. Therefore, we have from eq. ( p9|) 

Jd/2,9 = (d - 2)a d/2 ,g + qd(d- 2)o d / 2 -2, g -2, (63) 

^/2, g = (d- 2)ad /2>g + qd{d- 2)a' d/2 _ 2 q _ 2 + 4g(d - l)a d / 2 , 9 _2, (64) 

where the function afc jg , afc+i/ 2)(J and a' fc are some rea/ numerical constants 
listed in the Appendix. 

4.4 The total number of negative and zero modes 

Let us sum up the results of this section. 

q=l. In the case q — 1 there are negative modes of the operator A. There- 
fore, the effective action becomes complex. The imaginary part of the effec- 
tive action 'counts' the negative modes in the sense that 

\mT {x) = -\*N& i (65) 

where NyjH is the total 'number' of negative modes. Comparing ( |52|) and 
([65]) we find the number of negative modes to be 



JV YM 



f toUWntnH*, h = l). (66) 



a>0 



This causes, of course, the infra-red instability of the vacuum under small 
(quantum) fluctuations in the directions of the negative modes. 
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q=2, d=4. In the case q = 2, d = 4 there are no negative modes but there 
are zero modes. The effective action is real but it has a logarithmic infra-red 
divergency in the limit M — > 

r (1) (M)-= i<)log^ + O(l) ; (67) 

where Nyu is the 'number' of zero modes. Comparing eqs. ( |62"D and ( |67|) we 
find the total number of zero modes 



N (0) 

JV YM 



f dxJ2(^y 2 lGH 2 } (q = 2, d = 4) (68) 

The appearance of the logarithmic infra-red divergent term does not lead 
to serious physical problems. Whereas the zeta-function regularizes well the 
ultraviolet divergences, one has to renormalize additionally at the very end 
the coupling constants of the classical action in the infra-red region. This 
means that there appears an additional infra-red renormalization parameter 
A. In other words, after we have calculated Cym(0> M) there are no ultraviolet 
divergences any longer, but there are logarithmic infra-red ones. Thus one 
has to renormalize them as well. 



q=2, d>4 and q>3, d> 6. In the cases q = 2, d > 4 and q > 3, d > 6 there 
are no negative and zero modes. The effective action is real and the vacuum 
is stable at the quantum level, i.e., under small quantum disturbances. 



5 Analysis of the effective potential 

We define the total effective potential V(H) by 

r = s + r ( i) = fdxJ2 V ( H )- ( 69 ) 

Using the results of previous section with account of the classical Yang- 
Mills action calculated on the chromomagnetic background with equal mag- 
netic fields 



(70) 
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we find the general form of the effective potential to be 

V(H) = pi/ 2 + (a log * + b + c log A) , (71) 

where a, 6 and c are numerical constants depending only on the dimension 
of the spacetime d and the number of 'magnetic fields' q. 

The constants a, b and c are calculated in previous section. The form of 
the effective potential depends on the signs of these constants. The constant 
'a' is not equal to zero only for even d/2 > 2, i.e., in dimensions d — 4 (mod 4). 
The constant l V is always different from zero and is complex in the case q — 1. 
The constant c is not equal zero only in the case q — 2, d — 4. 

Since in the case q — 1 the effective potential becomes complex we do not 
study it in detail but concentrate our attention on the first nontrivial case 
q = 2 that can be investigated analytically to the end. 



5.1 q=2. (1=4 

In the case q = 2, d = 4 we have from previous section 

22 



a = 



3(4tt) 2 ' 
8 



3(47^ 



C(2) 



C + log2-log(27r) + l 



(72) 

(73) 
(74) 



Now one should absorb the infra-red divergency by renormalizing the 
coupling g. We define first 

M = Xe, (75) 

where A is arbitrary finite dimensionful parameter and e — * 0. Since V(H) 
cannot depend on the arbitrary parameters fi and A one assumes usually 
that the coupling constant depends on them in such a way that effectively 
the effective potential does not depend neither on /i, A nor on e, viz. 



(76) 
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where 

1 1 c A a A 

a(A) # 2 (/i, A) 4 A 2 4 /i 2 v 7 

^ = 7obH logE+6 - (78) 

" =a+c =-3(sp (79) 

and A is a constant scale of the 'magnetic field'. 

As a consequence we have two different renormalization group equations 
for g with respect to (j, and A 

^9 2 = M§), (80) 

X-^g 2 = f3m(g), (81) 

where 

M ) = -\? = -^\ (82) 

Ms) = ~f = (J^S 4 , (83) 



and an equation for a(A) = g 2 (A, A) 

, 9 ( 9 ^ d \-2 

A dA a= \^ + X dx) 9 



= Aot(«), (84) 

^=A=A 



Aot(a) = Ajv + Ar = -^« 2 = y^p« 2 (85) 
Thus we find important properties of the beta-functions 

Arv < 0, Ar > 0. (86) 

We see that according to the usual theory the effective coupling g(n, A) is 
asymptotically free in the high-energy limit fi — * oo, A being fixed. Moreover, 
it is asymptotically free in the limit A — > 0, \x fixed. 
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One should note that the behavior of the effective constant a (A) for 
A — > oo is different from the usual one. The contribution of the infra-red 
divergences (zero modes) changes the sign of the total beta function 



This means that the effective coupling a (A) is not asymptotically freel This 
effect occurs whenever the constant c is different from zero, i.e., only in the 
case q — 2, d — 4. In the standard case q — 1, d — 4 studied in the literature 
Ar = c = (there is no logarithmic infra-red divergency). Therefore, one has 
only one renormalization parameter fi and only one standard renormalization 
group equation fl8U| ) with /?uv < 0. This leads then automatically to the 
asymptotic freedom. 

The behavior of the effective potential is effected by the zero modes in a 
similar way. In our case, since (3 < 0, it is immediately seen (see Fig. 1) that 
although the vacuum H = is perturbatively stable the effective potential 
is unbounded from below and, hence, the perturbative vacuum is, in fact, 
metastable. 

One should note again that without the contribution of the zero modes 
one would have (3 = a > and, hence, the perturbative vacuum would 
be absolutely stable (see Fig. 2). These are zero modes that break down 
the stability of the vacuum at the nonperturbative level — the perturbative 
vacuum decays because of the zero modes. 



Consider now the case of odd dimension and odd d/2 > 3, i.e., d — 5, 7, 9, . . . 
and d = 6, 10, 14, . . .. Using the results of previous section we find that in 
both cases the constants 'a' and 'c' vanish 



Act > o. 



(87) 



5.2 q=2, d>4. 



5.2.1 d=2k+l, 4k+2 



a = c = 



(88) 



and the constant '6 



can be written in the form 



b 



sin(7r<i/4) 



{d-2f 



(Stt 2 )^ 



5 C(<*/2). 



(89) 



Ivan Avramidi: Stable chromomagnetic vacuum in higher dimensions 20 




Figure I: q=2, d=4, 5, 6, 7 (mod 8) 



d=8k+l, 8k+2, 8k+3. We see from that for d = 9, 10, 11 (mod 8), 
the coefficient h is positive and, therefore, the effective potential is positive 
and monotone (see Fig. 2). The vacuum H = gives the absolute minimum 
of the effective potential and is stable. 




Figure 2: q=2, d=9, 10, 11 (mod 8) 



d=8k+5, 8k+6, 8k+7. For d = 5, 6, 7 (mod 8), the constant b is negative. 
Qualitatively this case is similar to the case d — 4, q — 2. (See Fig. 1). Again 
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the effective potential is unbounded from below and the perturbative vacuum 
H = is metastable. 



5.2.2 d=4k 

Finally let us consider the case of even d/2 > 4, i.e., d — 8 (mod 4). The 
constant 'c' vanish 

c = (90) 

and the constants 'a' and '6' read 



n-( IV" 2<V2 "V- 2 ) 2 IP I 
a " (_1) (47r)2r(rf/2 + l) |Bd/21 ' 



(91) 



6 ^(li + z 1 2- c+log2 - log H- (92) 



d=8k+4. We see therefrom that for d — 8k + 4, i.e., d = 12 (mod 8), the 
constant a is negative. In this case the behavior of the effective potential 
repeats qualitatively that of the case d = 4. (Fig. 1) 



5.3 q=2. d=8k. 

The only case left is that of d = 8 (mod 8). This is the most interesting 
one. In all other cases there are no stable nontrivial, H ^ 0, vacuum states. 
Now the constant a is positive and this changes radically the behavior of the 
effective potential. 

First, it is bounded from below, i.e., there are some stable vacuums, at 
least the perturbative one, H = 0. It is not difficult to see that for some 
values of the parameters there is indeed some nontrivial local minimum that 
can also provide the absolute minimum of the effective potential. The form 
of the effective potential is governed by the coupling constant. (See Fig. 3). 
One can study the behavior of the effective potential in detail and find the 
following. 

There are three essentially different regions of the coupling constant 
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9\ ,9cv\t ,9 = 00 




Figure 3: q=2, d=8 (mod 8) 



I 

II 
III 



< g 2 < g\ 

2 „ 2 ^ 2 

9i< 9 < 9cr\t 

2 . 2 
5 >£crit 



(93) 



where 



2 8(d-4) 

0i = — - A exp 

da 



H, 



-(d-4)/2 



(94) 



where 



fl^ t = 2(d-4)^-(-- 4 )/ 2 , 

(X 



U 2 -i. 

H = fi e a . 



(95) 



(96) 



In the region I the effective potential is a monotone increasing function 
for H > and has only one minimum at H = — the perturbative vacuum. 
In region II there appear additional local minimum at some 



(97) 
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where h m - m is the solution of the equation 



(98) 



with 




(99) 



The values of if m i n lie in the region 



Hi < H mm < ifcrit, 



(100) 



where 



Hi = H exp 



4(d-2) 
d{d - 4) 



(101) 



2 



) 



(102) 



f/crit = H exp 



d-4 



This local minimum indicates that there is an additional metastable vac- 
uum state with the energy V m \ n larger than that of the perturbative level, 
Kim > 0. Such metastable state must decay in the usual perturbative vac- 
uum state H = 0. 

At the critical value of the coupling constant g = g crit the effective poten- 
tial has the minimum at H min = H crit and the energy of this additional state 
becomes equal to that of the perturbative vacuum, V m i n = 0. This means 
that both states, H = and H = H CT n, are physically equivalent. The true 
vacuum state in this case is a mixture of the perturbative vacuum and the 
nonperturbative one. In other words, at the critical value of the coupling 
constant a phase transition occurs: a new stable vacuum state appears and 
the system decays in it. 

In region III the perturbative vacuum is either metastable or absolutely 
unstable. The minimum of the effective potential is provided by the new 
nonperturbative stable vacuum state H = H vajC . The values of H vac range in 



-f^crit < H vac < if 2, 



(103) 



where 



H 2 = H e d . 



(104) 
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It is interesting to note that the values of the nontrivial vacuum fields are 
bounded from above by H 2 — there are no vacuum states with if vac > H 2 . 
Besides, as the coupling constant g increases, the vacuum energy decreases 
and reaches its absolutely minimal value 

min Kac = -^# d/2 . (105) 

0< 9 <«d de 

6 Concluding remarks 

In this paper we continued the study of the low-energy behavior of the non- 
Abelian gauge theories initiated in [[|. We applied the covariant algebraic 
method for the calculation of the heat kernel of Laplace type operators de- 



veloped in [TR [16| to the case of covariantly constant background fields, 



with q 'magnetic fields' of equal amplitudes, (see eqs. ([18]), (|T9|), (|25|) ). The 
heat kernels, zeta-function and the one-loop effective action are calculated 
explicitly in terms of the function b s ^ q (z) that is studied in detail in the 
Appendix. The cases q = 1 and q = 2 for any spacetime dimension d are 
studied analytically to the end — the zeta-functions and the effective action 
are calculated exactly. 

We confirmed the old result that in the case of only one 'magnetic field', 
q = 1, the vacuum is essentially unstable, i.e., unstable under small (quan- 
tum) disturbances. In other words, there are negative modes that destroy 
the vacuum immediately — there are no stable vacuum states with only one 
magnetic field. 

In the case of two 'magnetic fields' of equal amplitudes the effective po- 
tential is studied in detail. We distinguish two different cases d = 4 and 
d > 5. For d = 4 there are zero modes of the gauge fields. These zero 
modes do contribute to the renormalization group beta functions and rad- 
ically change the asymptotic behavior of the effective coupling constant - 
it is not asymptotically free any longer. The perturbative vacuum H = 
is stable under small (quantum) disturbances but is unstable under classical 
variations of the background field. 

For d > 5,6, 7, 12 (mod 8) the vacuum is again metastable (Fig. 1). For 
d = 9, 10, 11 (mod 8) the perturbative vacuum H = is absolutely stable 
(Fig. 2). For d = 8 (mod 8) and sufficiently large coupling constant there ap- 
pears some nonperturbative stable vacuum state H = if vac ^ with negative 



Ivan Avramidi: Stable chromomagnetic vacuum in higher dimensions 25 



vacuum energy V vg , c < (Fig.3). 

One should note that our analysis is applicable directly to Euclidean field 
theory. In case of spacetime of Lorentzian signature there is a more strong 
restriction d > 2q + 1 — there must be enough space to put q magnetic 
fields. That is why our results are not applicable to the case of Lorentzian 
spacetimes of dimension d = 2q ||, e.g. for q — 2, d — 4. 
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Appendix. The function b SA (z) 

In this appendix we study the function b s ^ q (z) defined by 

1 00 -tz 
v ' 

It is not difficult to show that for Re z > —q, the integral ( |106| ) converges 
in the region Res < and defines an entire function of s by analytical 
continuation. In particular, in the region Res < N, N being some natural 
number, it is defined by 



bs,q{z) 



-1 



r(-s + iv) 



dtt 



-s-l+JV 



-tz 



( 1 



\ sinh t 



(107) 



For non-positive integer s ^ 0,1,2, ... the function b s ^ q (z) as a function 
of z has a branching point z = —q. It is analytic function of z in the complex 
plane cut along the real axis from z = —q to — oo. 
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in 



At positive integer points s = 0,1,2,..., exactly in the same manner as 
we find that the values of the entire function b Stq (z) are determined 



by the Taylor expansion of the integrand 



dt 



-tz 



t \ q 



sinht 



(A; = 0,1, 2, 



(108) 



t=o 



For half- integer s — k + 1/2, by choosing N = k + 1 in ( p. 07]) we get 



k+l 



-* ( t y 

\ sinh t 



(k — 0, 1, 2, 



(109) 

By differentiating eq. ( p. 07] ) and choosing iV = k + 1 we also obtain the 
derivative 

d 

b s, q {z) = -Q- S h ^i( z ) 



at integer points s = k 

b' k , q {z) = -Cb k , q (z) - j dt log tl-- 



k+l 



-tz 



V sinh t 



(110) 



For nonpositive integer q the function b s ^ q (z) is expressed in elementary 
functions 



r -s) 



0<n<-q 



-l) n (z + q + 2n) s - q 
n\ r(— q — n + 1) 



(111) 



In particular, 



(q = 0,-1,-2,...). 



K-i(z) 



2(s + 1) 



,2 + l) s+i - (z - 1 



,s+l 



(112) 
(113) 



etc. 
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For q > 1 one can expand 1/ sinh 9 t in the powers of exponents to get an 
infinite series representation 



b s Jz) = 2«(-l) 



r(s + i) 



T(s-q + l)T( q ) „> 
Using the function (see [^3], p. 27]) 

oo 

1 



E^rr^(* + </ + -M*-" (Hi) 



n! 



$(A,p,u) 



r(p) 



1 - Ae- 



E 



A n 



(115) 



one has finally 
&.,„(*) = (-1)«2* 



r(s + i^ 



5 



r(s-g + i)r(?) v^a / 



A ?-1 $(A,g-s. 



A=l 

(116) 

Some properties of the function $(A,p, u) as well as its connection with the 
Riemann and Hurwitz zeta-functions are described in | 23| . 

Let us study now the boundary value of the function b StQ (z) and its deriva- 
tives as z — > — is 



clef 



lim b sn (z) 

z—fO—ie 



(117) 



The values of the function b Stq (z) at positive integer points given by ( |10§| ) 
are polynomials in z 



where 



Ofc, 



h,q(z) 



dt 



E 

0<n<fc 



t 



z k n a. 



n,qi 



sinht 



(£ = 0,1,2,...) 



;il8) 



(119) 



t=o 



Taking the limit z — > and observing that the Taylor series of the function 
(t/ sinht) 9 is a power series in (— t 2 ) with positive coefficients we find 



a 2 k+i, q = (A; = 0,1,2, 



(120) 
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and the constants a2k,q have the property 



0>2k, 



■l) k \a 2k , 



;i2i) 



The values of function a S}Q at half integer points s = k + 1/2 and its 
derivative a' s at integer points s = k depend crucially on q. For nonpositive 
q the function b s ^ q (z) is not analytic at the point z — 0. Therefore, one has 
to take into account the infinitesimal imaginary part in z — > — is. 

We have for q — — 1 



1 



and for g = 



1 

~2(fc + l) : 



a fc+l/2,0 — Ofc+1,0 



2A; + 3 



1 + <(-!) 



fe+i 



+ 



-l) fc+1 i7T 



2(fc+ 1)' 



(122) 
(123) 



(fc = 0,l,2,...) 



For = 
singular 



a' k q is not well defined because the derivative b' s 



log z + o(z). 



(124) 
z) at s = is 
(125) 



As we already mentioned above b Sjq (z) is an entire function of s for Re z > 
—q. This means that for q > 1 the function b Syq (z) also remains analytic in 
the limit z — > 0, i.e., a Sj(? is entire function of s for g > 1. Therefore, for 
q > 1 we can just put z = in eqs. (|109|) and ( |1 10|) to obtain 



and 




(126) 



(127) 



(A; = 0,1,2,...). 
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We see that ak+i/2, q and a' k q are some real numerical constants. 
There are two important particular cases when the integral 



T(-s) 



dtt 



-s-l+q 



1 



sinh 9 1 



;i28) 



can be calculated analytically for q > 1 too, namely, q = 1 and q = 2 
Substituting z = in eq. (|114| ) we obtain for q — 1 

a Stl = -2s(l-2 s ~ l )C(-s + l) 



and for q = 2 



a s , 2 = 2 s s(s-l)C(-s + l). 



Using the property [23 



ai-s) = 2(2n)- s cos[-s)T(s)as) 



7T 



(129) 
(130) 

(131) 



we obtain the values of the functions a Sj i and a s _ 2 and their derivatives at 
half-integer points and integer points 



a k+l/2, 1 



\[(*+l)/2] 



2 3/2( 2 fc-i/ 2 _ i)(2 7r )- fc -V2 F ( A . + 3/2)C(A; + 3/2) (132) 



l 2fc+l, 1 



l) fc + 1 (2 2fc - l)(2fc + l)!(27r)- a *C(2A; + 1 



-2fc, 



(133) 



4 2fc, 1 



C(2k) 
C(2k) 



+ $>(2k + l) + 



)2fe-l 



2 2fc-l _ ! 



log2-log(2vr) 



(134) 



a fc+ i/ 2 ,2 = (-l) [(fc+1)/2] (2A; - I^tt)-*" 1 / 2 !^ + 3/2)C(* + 3/2) 



l 2k+l,2 



-l) k+1 2 2k 2k(2k + l)!(2vr)- 2fc C(2A; + 1) 



(135) 
(136) 



*2fc, 2 



^-l) fc 2 2fc+1 (2A; - l)(2fc)!(27r)- 2fc C(2A:) 



x 



^ + *(2* + l) + lo g 2-lo g (2.) + ^ T 



(137) 
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